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Abstract 

Bose-Einstein condensation of charged scalar and vector particles may actually oc- 
cur in presence of a constant homogeneous magnetic field, but there is no critical 
temperature at which condensation starts. The condensate is described by the sta- 
tistical distribution. The Meissner effect is possible in the scalar, but not in the 
vector field case, which exhibits a ferromagnetic behavior. 
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It was pointed out long ago by Schafroth |]T| that for a non-relativistic boson gas, Bose- 
Einstein condensation (BEC) does not take place in the presence of a constant magnetic 
field. The problem was studied afterwards by May , who investigated the condensation 
in the non-relativistic case for an arbitrary dimension d and showed that it occurs for 
d > 5, and later by Daicic et a/ 0, who extended the considerations made by 
to the relativistic high temperature case and found new magnetization properties. More 
recently Toms has proved that BEC in presence of a constant magnetic field does not 
occur in any number of spatial dimensions, and Elmfors et al consider that a in the "id 
case, although a true condensate is not formed, the Landau ground state can acommodate 
a large charge density. The last paper discuss also the magnetization of the relativistic 
scalar gas and shows that the Meissner effect occurs in the low temperature case, in 
analogy to the non-relativistic case studied by All these papers are characterized 
by very careful calculations and their discrepancies, when arised, are in general related 
to subtle points. One essential argument in all of them is that condensation cannot 
occur since the condition for condensation, taken from the symmetric (zero field) case 
as /i = V Mc^ ± eBfic for the relativistic case {fx' = fx — Mc^ ± eBfi/mc = in the 
nonrelativistic limit), when applied to the case in which there is an external magnetic 
field, leads to a divergent behaviour of the density in terms of /i'. 

There are two different ideas which usually are considered to be the same, concerning 
what is to be understood as BEC: 1) The existence a finite fraction of the total particle 
density in the ground state at some temperature T > 0. 2) The existence of a critical 
temperature Tc > for which /x' = 0, such that for T < Tc some significant amount of 
particles start to condense in the ground state. 

The present author's point of view is that i) the condition for condensation in the 
magnetic field case corresponds to 1), i.e. it cannot be extrapolated as the same one 
of the standard zero field case, (which is 2)), since in the magnetic field case there are 
different physical conditions: explicit spatial symmetry breaking, discrete Landau states; 
ii) since the chemical potential is not an independent thermodynamical variable, but 
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depends on the density, temperature and magnetic field, there is no divergence problem at 
all, iii) as different from the zero field case, in presence of the magnetic field the expresion 
for the density contains the ground state contribution. 

Let us remind the origin of the critical quantities /ic, Tc in the standard theory of 
BEC (with zero magnetic field). The chemical potential fi' = f{N,T) is a decreasing 
function of temperature at fixed density A^, and for fi' = one gets an equation defining 
Tc = fc{N). For temperatures T < as /i' = 0, the expression for the density gives values 
N'{T) < N, and the difference — A^' = A^^o is interpreted as the density of particles 
in the condensate. This is to be expected since the density of states, proportional to 
Anp"^, cancels the infrared divergence of the Bose-Einstein distribution n{p,T) for /i = 0. 
The expression for A^ contain only the particles in excited states. Thus, for /i' 0, 
taking pq ~ a/— 2Myu', the amount of particles in a neighbourhood of the ground state 
of amplitude Pq is A^o = 47r Jq" p'^n{p,T)dp ~ MTpo/h^. In the magnetic field case, the 
infrared divergence of the Bose- Einstein distribution resulting from taking yu' = is not 
cancelled by the density of states in momentum space at T 7^ 0. In other words, for a 
given constant density, we are not allowed to put fi' = keeping T 7^ 0. Now the amount 
of particles in a small neighbourhood of the ground state of amplitude po, as we shall 
see below, is A^o ~ /q" n{p2,,T)dp3 ~ eBM^^'^T/ y/^JF, which is large for small p.' and 
leads to conclude that in the magnetic field case the expression for the density accounts 
for the contribution of the population in all quantum states. 

Thus, in the magnetic field case the problem of finding a macroscopic ground state 
density is not conditioned to have p' = 0. 

We shall investigate in more detail the problem of BEC in the physical d = 3 space, 
at large densities and strong magnetic fields. Under these conditions BEC occurs, but 
without having a definite critical temperature Tc > 0, i.e. the phase transition is diffuse. 
The concept of diffuse phase transitions, as those not having a definite critical tempera- 
ture, but an interval of T, was introduced long ago by Smolenski and Isupov in their 
investigation of the phase transition which occur in certain ferroelectric materials. 
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By assumming as in a constant microscopic magnetic field B along the ps axis 
(the external field is H'^^^ = B — AnAA^B), where Ai{B) is the magnetization), the 
thermodynamic potential for a gas of scalar particles placed in it is 



eB 



= / dp, ln(l - e-(^-'^)^)(l - e-(^'+^)^) + /3e, (1) 

A-K^n Cfj J -oo L J 

where = \Jp1c^ + M^c^ + 2eBTic{n + |), the last term in (1) accounts for the vacuum 
energy and /i is the chemical potential. We shall use (3 = 1/T, where T = kt m what 
follows, k being the Boltzmann constant and t the absolute temperature in Kelvins. 
For a vector field the one-loop thermodynamic potential is 

eB 



a 



dp, ln(l - e-(^°'-^)^)(l - e~(^«'+^)^) + (3e^ 



+ 



eB 



^ /3„ / ln(l - e-^^'-'^)^) (1 - e-(^^+^)^) + (5e, 



(2) 



where /3„ = 3 - 5on, eog = V^ic^ + ^^c^ - e5/ic, = ^plc^ + M2c4 + 2eBhc{n + i). 

The mean density of particles minus antiparticles (average charge divided by e) is 
given by N^^^ = -dQs,v/dfi. 

We have explicitly 
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where rip = [exp{eq ^ fi)l3 — 1] ^. For the vector field case we have, 
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where n^p = [exp{eoq =f - eog = -y/PsC^ + ^^c^ - e5^c, and /?on = 3 - Son- 

For T ^ 0, /i ^ M±c^ (we named M± = J M"^ ± eBh/c^), the population in Landau 



quantum states other than n = vanishes (this is shown explicitly in the Appendix) 
and the density for the n = state is infrared divergent. We expect then most of the 



population to be in the ground state, since for small temperatures n^p is vanishing small 
and bell-shaped curve with its maximum at ps = 0. We will proceed as in H| 

and call po{'^ y/—2M fx') some characteristic momentum. We have then, by assumming 
—fi' <^ T, for the density in a small neighbourhood of ps = 0, 



eBT fpo dps 



2n^h^c Jo ^p2^2 + ]vf2c4 ± eBhc - /i 
eBT fPo dps 



eBT 2M± 

(5) 
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where Ng^v = Nos,v + SNg^v and 6Ns^v is the density in the interval [po, oo]. Actually 6Ns^v 
is negligibly small and N^^^ ~ A^s,«- 

The expression (|^) was obtained in p[| and more recently in ^j. Formally it indicates 
a divergent behaviour of A^os,i; for A^' 0, but we must be careful in doing that interpre- 
tation. Actually it means nothing more that in a neighbourhood of the ground state, at 
constant temperature, A^os.i) ~ and for constant Nqs^v, yu' ~ T'^- This comes from 

the fact that /i' is not, and cannot be taken, as an independent thermodynamic variable 
(see i.e. [0, chapter 3), but depends on the quantities N^^,T,B and it is enough that 
fi' ~ T'^f{Nos,v) for T (/ being a finite function) to avoid the divergence. Thus (H) 
simply means that 

, _ e'B^T'M^ 

We observe that yu' is a decreasing function of T and vanishes for T = 0, where the 
"critical" condition /i = M±c^ is reached. As shown below, in that limit the Bose-Einstein 
distribution degenerates in a Dirac 5 function, which means to have all the system in the 
ground state p^ = 0. To see this, we shall rewrite the momentum density of particles 
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around the ground state pa = 0, n = approximately as 

T 



2Af± + 87r2a*c27V^^„ 



where 



where pr = \f2jiM±T is the thermal momentum, = hcX/eB the elementary volume 
cell, A = h/pT being the De Broglie thermal wavelength. We have approximated the Bose- 
Einstein distribution by one proportional to a Cauchy distribution, having its maximum 
~ 7^^ for p3 = 0. We have that 7 ^ for T ^ 0, but for small fixed T, 7 also decreases 
as v^Nos,v increases. We remind that in the zero field case, the condensation condition 
demands A^A^ > 2.612. 
One can write then 

1 eB n . . 

2^0s,v = ^^2f2^ no{p3)dp3. (8) 

Thus, approximately one half of the total density is concentrated in the narrow strip 
of width 27 around the Ps = momentum. It results that for densities and magnetic 
fields large enough, if we choose an arbitrary small neighbourhood of the ground state, of 
momentum width 2p3o, one can always find a temperature T > small enough such that 
7 ^ Pso and (|^) and (H) are satisfied. The condensate appears and it is described by the 
statistical distribution. 
We have also 



and obviously 



limno{p3) = A7f^h^c^6{p3). (9) 
7^0 eB 



roc 

Ns^v = lim^^Q o / no{ps)dp3. (10) 
A-K^h C J-00 
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Thus, if T = 0, all the density Ns^v lies in the condensate, as occur in the zero field case, 
but here the total density is described explicitly by the integral in momentum space. 

Obviously, one cannot fix any (small) value for 7 from which the distribution starts 
to have a manifest S{p3) behavior; and there is no critical temperature for condensation 
to start, which we interpret as a diffuse phase transition. But when 7 and 1/A decrease 
enough (we define A = uh/T, and uj = cB/Mc. If A ^ 1, the system is confined to 
the n = Landau state, see below), the conditions for condensation mentioned above are 
satisfied. 

In the non-relativistic case, the above results can be derived even from Schafroth's 
formulae. We can write the non-relativistic limit for A^^ (which is a very good approxi- 
mation for the relativistic case if Mc^ ^ T, since the main contribution to the integrals 
comes from values of < Mc ) , as 



eB nMT ~ e"^!'"^ 1 



2n^h^c^ 2 ^1 ml/2 i _ ^-2u.nmi3 



eB InMT g e"'^!'"'^ 



2rf^V 2 ml/2 



1 + 



^—2uihm/3 



m=l 

where the unity in angular brackets accounts for the Landau n = state and the second 
term for all the set of excited states. We have used /ii = Mc^ + uTi — /i and u = eB /Mc. 
Now, if the parameter A = uh/T ^ 1, and if /ii — > 0, even at large temperatures, we 
can neglect the second term in angular brackets in (0), and approximate the resulting 
sum by an integral. This is done by introducing the continuous variable x = Tm and by 
writing TYl,'^=i Iq^ dx. After integrating, we obtain back (5), with /ii, M in place of 

Neutron stars, where strong magnetic fields and very high densities are assumed to 
exist, may provide conditions for the ocurrence of BEG. It has been conjectured that 
superfiuid and superconductive effects are produced (see i.e. [|^, and references 
therein). 

As a simplified model of the gas of kaons in a neutron star, we take A^^ ~ 10'*'^ cm~^. 
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t ~ 10®°K, and local magnetic fields B ~ 10^^ G (fields of order 10^^ G have been estimated 
inside hadrons the conditions for BEG can be also satisfied, if the medium provide 
also screening mechanisms for the very large electric fields which also arise. In this case 
7 = 10^^°. If we suppose the star of dimensions ~ lO^cm, the discrete momentum states 
would be spaced by an amount 6p ~ lO^^^gcm/s. This means to have one half of the 
total density distributed in these 2'-f/6p = 10^ quantum states. The ground state density, 
with strictly zero momentum p^, can thus be estimated as 10'^'' cm~^. We observe in this 
case the dimensionless phase space density Nv^ ~ 10^^. 

We turn now to the magnetization problem. From (|l]), we get that for T ^ 0, the 
magnetization of the scalar field is 



E/ rfp3^— + (12) 





In the condensation limit we get 

tl- 

We see that the magnetization is opposed to the external field and the system behaves as 
a perfect dia-ferromagnetic (this was first pointed out by Schafroth 0]). 

If we take the B{> T) limit of (|1^), one can discuss the critical conditions for 
the arising of the Meissner effect in the low temperature relativistic case. The condition 
(obtained for T 7^ in ref. 0), is H^^^ = —A4s{0). In particular, we agree with that 
in the relativistic case the Meissner effect occurs in analogy with the non-relativistic one, 
and that it is not neccessarily connected with high temperature pair creation proccesses, 
as suggested in [Q. 

For the vector field case, the magnetization in the condensation limit is positive since 
all the system is in the Landau n = state, (see eq. (2)), and 

. , eN„h , , 



and we have that the condensate of vector particles behave as a true ferromagnet. In 
particular, we can write = H^^^ = B — A4{B) as the condition for spontaneous magne- 
tization to occur. 

In concluding, it is important to remind that for the first time true BEC has been 
recently observed [|I^ in evaporated ^^Rb atoms confined in a magnetic trap creating 
an ellipsoidal potential. The system is far from being an isotropic noninteracting gas, 
and although it also differs from our present model, it has, however, some analogies. In 
particular, the anisotropy of the confining potential leads to a larger velocity spread in 
the axial direction as compared to the radial one, which is the analog of the ps spread of 
the n(p3) distribution in the magnetic field case, confined radially to the n = Landau 
quantum state. 

The author thanks R. Baquero, J. Hirsch, K. Kirsten, D. B. Lichtenberg, O. Perez- 
Martinez and D. Quesada for comments. He is especially indebted to A. E. Shabad for a 
discussion long ago, from which some of the basic ideas presented in this paper arised. 



We will present a demonstration concerning the vanishing of the density for excited Lan- 
dau states if A ^ 1, Mj3 ^ 1. We use in this section units h = c = 1, and consider only 
the scalar field case. Let us call the density corresponding to excited states. We can 
write, by using the K2 Bessel function representation in terms of dimensionless quantities 
M_i_ = M_|_/?, M+ = and fl = fif3, after summing over Landau quantum numbers 

from 1 to cxD, 



We will establish an upper bound to this quantity. Let us introduce x = M^t and cut 



Appendix 




(15) 



the integral in two parts by the point Xq = 6/ (M_|_ + A), where ^ ^ 1. Let us call also ^ 



to some point in between xq and 00. We have 



8 



Both series can be made to converge to an arbitrary small number as 6' — > and 
1/A ^ 0, provided that M > 1. 
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